model:
sets:
    area/1..7/: POP, X, Y;
endsets
data:
    POP = 38,71,51,32,24,49,66 ;
enddata
max = @sum(area(i): POP(i) * Y(i));
Y(1) <= X(1) + X(2) + X(3)+X(7);
Y(2) <= X(1) + X(2) + X(3);
Y(3) <= X(1) + X(2) + X(3)+X(4)+X(6)+X(7);
Y(4) <= X(3) + X(4) + X(5)+X(6);
Y(5) <= X(5) + X(6) + X(4);
Y(6) <= X(3)+X(4) + X(5) + X(6)+X(7);
Y(7) <= X(1) + X(3) + X(6) + X(7);
@sum(area(i): X(i)) = 2;
@for(area(i): @bin(X(i)));
@for(area(i): @bin(Y(i)));
end




model:
sets:
g /1..7/:a;
coo(g,g):x;
endsets
max = 109*x(1,2)+104*x(1,3)+89*x(1,4
)+122*x(2,4)+117*x(3,4)+115*x(3,5)+100*x(4,5)
+83*x(4,6)+81*x(5,6)+73*x(5,7)+56*x(6,7);
@for(coo(i,j):@bin(x(i,j)));
x(1,2)+x(1,3)+x(1,4)+x(2,4)+x(3,4)+x(3,5)+
x(4,5)+x(4,6)+x(5,6)+x(5,7)+x(6,7)=2;
x(1,2)+x(1,3)+x(1,4)<=1;
x(1,2)+x(2,4)<=1;
x(1,3)+x(3,4)+x(3,5)<=1;
x(1,4)+x(2,4)+x(3,4)+x(4,5)+x(4,6)<=1;
x(3,5)+x(4,5)+x(5,6)+x(5,7)<=1;
x(4,6)+x(5,6)+x(6,7)<=1;
x(5,7)+x(6,7)<=1;
end

